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摘要 利用 中 国 科 学 院 新 疆 天 文 台南 山 观测 站 26 m 射 电 望远镜 , 在 中 心 频率 1556 MHz, 对 Crab 脉 冲 星 (PSR 
B0531 上 +21) 进 行 了 长 达 12.6 h 的 连续 观测 , 观测 带宽 为 512 MHz, 时 间 分 辨 率 为 32 us, 研究 了 巨 脉冲 辐射 的 等 待 
时 间 分 布 特征 . 观测 共 探测 到 2097 个 信 噪 比 大 于 10 的 巨 脉冲 , 对 应 的 流量 密度 大 于 100 Jy. 巨 脉冲 的 爆发 率 表 现 
FARE UTA ME, 在 较 短 的 时 间 内 具有 较 高 的 爆发 率 , 在 相对 长 的 宁静 期 内 巨 脉 冲 的 爆发 率 较 低 , 尤其 是 中 间 肪 
冲 相位 内 的 巨 脉冲 爆发 . 相 邻 两 个 巨 脉冲 的 等 待 时 间 分 布 表 现 为 寡 律 分 布 特征 , 可 以 用 一 个 非 稳 态 的 泊 松 过 程 
进行 模拟 , 这 表明 巨 脉冲 的 爆发 是 一 种 独立 的 随机 事件 . 此 外 , 主 脉冲 和 中 间 脉 冲 相 位 上 的 巨 脉 冲 有 具有 不 同 的 等 
PER E, 征 , 这 意味 着 脉冲 星 不 同 磁极 的 巨 脉冲 辐射 机 制 可 能 是 不 同 的 . 这 些 观测 结果 对 于 理解 脉冲 星 的 
射电 辐射 机 制 要 意义 . 

关键 词 ”脉冲 星 : 个 别 : Crab, 脉冲 星 : 巨 脉 冲 , 方法 : 观测 , 方法 : 数据 分 析 
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1 引言 后 被 称 为 巨 脉冲 (Giant pulse, GP)". GP 出 现在 

脉冲 星 是 快速 旋转 的 、 高 度 磁化 的 中 子 星 , 在 所 有 射电 频率 上 , 并 只 存在 于 主 脉 冲 (Main pulse, 
它 的 磁极 处 产生 电磁 辐射 束 . 当 脉冲 星 旋转 时 , 地 MP) 和 中 间 脉 冲 (Interpulse, IP) 相 位 内 的 辐射 窗 
球 上 的 射电 望远镜 就 会 接收 到 周期 性 的 单个 脉 溃 口中 , 与 MP 相 比 卫 窗 口中 出 现 的 频率 大 约 低 一 个 
辐射 . 由 于 单 脉冲 辐射 非常 微弱 , 因此 大 部 分 脉 ” BCE ARP. GP 是 宽带 辐射 、 持 续 时 间 短 、 时 间 尺 
冲 星 的 单 脉冲 很 难 被 探测 到 . Staelin 等 人 在 1968 年 ”” 度 只 有 几 微 秒 , 结构 时 间 尺 度 可 达 10 nsl. 脉冲 
首次 观测 到 来 自 于 Crab 脉 冲 星 极 强 的 单 脉冲 辐射 ， 流量 超过 平均 脉冲 的 10 倍 以 上 , 能 量 遵循 守 律 分 
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MA, 关于 GP 起 源 的 研究 已 经 投入 了 大 量 的 工作 ， 
例如 Mikhailovskii 等 回 认 为 GP 是 非 线性 孤立 子 通 
过 脉冲 星 磁 层 等 离子 体 传播 所 产生 的 . Weathe- 
ral 认为 强 辐射 是 | 含有 高 激发 静电 灌流 的 相对 
论 等 离子 体 产生 的 . Petrova” tE HGP Æ HEME 
普 顿 散射 导致 无 线 电 发 射频 率 的 重新 分 布 产 生 的 . 
Istominl 引 认为 GP 是 由 脉冲 星光 速 圆柱 区 域 的 磁力 
线 的 磁 重 联 放 1 


1 年 


用 于 GP 加、 模式 变换 PI、 脉冲 消 零 史 和 快速 射 
RS EI 究 . Lundgren 等 26] 发 现 Crab 脉 冲 星 GP 
的 到 达 时 间 服 从 泊 松 统计 ， 与 GP 相 邻 的 脉冲 周 
期 与 任何 其 他 脉冲 周期 具有 相同 的 GP 辐 射 概率 . 
Karuppusamy 等 站 得 出 Crab 脉 冲 星 中 GP 的 到 达 时 
间 也 呈 泊 松 分 布 , 这 进一步 表明 GP 活动 的 存在 既 
不 抑制 也 不 增强 辐射 . 然而 , 在 实际 观测 结束 时 ， 
由 于 灵敏 度 降低 , GP 的 到 达 时 间 最 终 呈 现 出 长 尾 


EP ÆJ. Lyutikov 卜 提出 GP 是 1 
在 靠近 光速 圆柱 的 封闭 磁力 线 上 通过 反常 回旋 共 
振 产 生 的 . 由 于 GP 具有 极 短 的 爆发 时 标 , Hankins 
等 四 认为 GP 起 源 于 非 热 的 相干 辐射 过 程 . 以 上 的 
模型 预测 都 无 法 完全 解释 观测 , 所 以 GP 的 辐射 机 
制 没有 定论 . 
到 内 外 研究 团 


组 已 经 对 Crab 脉 冲 星 的 GP 开展 
了 大 量 的 观测 研究 Ll. 在 较 低 观测 频率 , GP 在 传 
番 的 过 程 中 会 受到 星际 介质 中 等 离子 体 的 影响 , 产 
E 色散 和 散射 效应 , 脉冲 出 现 明 显 展 宽 , 并 且 GP 的 
# 量 表现 出 寡 律 分 布 特征 112. Sallmen 等 1 引 在 射 
电波 段 对 Crab 脉 冲 星 的 GP 进 行 了 多 波段 观测 , 探 
则 到 了 29 次 GP 爆发 现象 ,得 到 了 GP 的 谱 指 数 范 
围 为 -2.2 到 一 4.9. Hankins 等 四 在 5.5 GHz 频率 上 对 
Crab 脉 冲 星 的 GP 进 行 了 观测 , 发 现 GP 包 含 多 个 持 
续 时 间 在 纳 秒 量 级 的 子 脉冲 辐射 结构 , 这 与 Wea- 
theralltg 提 出 的 等 离子 体 扰动 模型 相符 合 . GP 的 
能 量 和 脉冲 宽度 具有 相关 性 , 持续 时 间 越 短 的 脉冲 
对 应 的 爆发 能 量 越 强 ,1 引 ，Cordes 等 中 对 Crab 脉 
冲 星 的 GP 进行 了 多 波段 的 观测 , 发 现 GP 仅 在 MP 
HIP 的 相位 处 出 现 , 在 其 他 轮廓 成 分 相位 处 未 探测 
1GP 爆 发 , 在 射电 中 高 频率 上 , IP 相 位 处 的 GP 爆 
发 逐渐 占 主导 . Karuppusamy 等 外 在 MP 和 IP 相 位 
处 同时 探测 到 了 GP 辐射 , 并 且 双 GP 的 性 质 与 单 
个 出 现 的 GP 之 间 没 有 区 别 . 通过 高 分 辨 率 观测 ， 
Jessner 等 Li" 探测 到 了 更 加 微弱 的 爆发 现象 , 发 现 
位 于 MP 和 IP 相 位 处 的 GP 在 轮廓 形状 和 偏振 特性 
方面 存在 巨大 的 差异 , 这 表明 它们 在 磁 层 中 的 传播 
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分 布 ( 威 布 尔 分 布 ). 据 报道 , 毫秒 脉冲 星 B1957 二 20 
的 GP 可 以 用 威 布尔 分 布 更 好 地 描述 , 并 且 可 能 
辐射 模式 的 变化 有 关 B271. PSR J1326-6700 模 式 变 
化 时 间 尺 度 的 内 在 分 布 也 可 以 用 威 布尔 分 布 很 好 
地 描述 , 它 在 模式 变换 中 具有 一 定 程度 的 记忆 效 
应 B3. 巨 微 脉冲 的 WTD 也 表明 Vela 脉 冲 星 具 有 
定 程度 的 记忆 效应 Pa. 以 往 的 观测 虽然 探测 到 了 
大 量 的 GP 样本 , 但 是 由 于 观测 不 连续 , 所 以 无 法 
统计 GP 的 等 待 时 间 . 本 文 需要 对 巨 脉冲 的 等 待 时 
间 进 行 统计 分 析 , 所 以 需要 进行 长 时 间 的 连续 观 
M. 通过 计算 相 邻 两 个 GP 之 间 的 时 间 间 隔 ( 等 待 时 
间 , WT), 可 以 获得 可 靠 的 统计 数据 , 这 可 以 提供 关 
于 GP 是 作为 独立 事件 发 生还 是 由 其 他 事件 连接 / 
触发 的 信息 , 并 进一步 为 GP 辐射 模型 提供 观测 依 
据 . 


本 文 分 析 了 Crab 脉 冲 星 1556 MHz 频 率 GP 的 
WTD, 并 用 非 平 稳 泊 松 过 程 的 统计 模型 进行 拟 合 . 
文中 第 2 节 介 绍 了 观测 和 数据 处 理 过 程 ; 第 3 节 和 
， 对 观测 结果 进行 了 详细 分 析 和 讨论 ; 第 5 

是 论文 研究 的 总 结 . 


2 ”观测 和 数据 处 理 

在 2020 年 10 月 17 日 15:58:09 UT, 用 南山 26 m 射 
有 望远镜 (Nanshan 26-m Radio Telescope, ae 
对 Crab 脉 冲 星 的 单 脉 冲 进行 了 长 达 12.6 h 的 连续 
测 . NSRT 接 收 机 由 两 个 线 偏 振 
其 中 心 频率 为 1556 MHz、 带 宽 512 MHz, 其 中 约 


= 


路 径 也 是 不 同 的 . 

等 待 时 间 分 布 (Waiting time distribution， 
WTD) 在 地 球 物理 和 天 体 物理 学 中 有 着 广泛 的 应 
gas 20. 在 脉冲 星 物理 学 中 , 人 们 一 致 认为 WTD 
DWH ae IS eo ti, WTD 在 过 去 十 年 中 一 直 被 
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1/3 带 宽 受 到 射频 干扰 (Radio Frequency Interfer- 
ence, RFI). 接收 到 的 信号 首先 通过 降 频 变换 和 放 
大 , 接着 被 传递 到 基于 ROACH2 (Reconfigurable 
Open Architecture Computing Hardware Genera- 


tion 2) 的 数字 多 相 滤 波 器 组 系统 后 端 , 生成 512 个 
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SUE IS). 然后 系统 以 32 hs 的 有 效 采 样 时 间 、 
8-bit 的 精度 对 数据 流 进行 采样 , 为 了 方便 后 续 线 下 
处 理 , 采样 数据 以 PSRFITS (Pulsar Flexible Image 
Transport System) 格 式 保存 9. 

在 线 下 数据 处 理 中 , 首先 对 频率 通道 逐个 进 
行 统计 分 析 , 搜索 显著 偏离 正 态 分 布 的 通道 , 然后 
去 除 相 对 应 的 时 间 序 列 和 频率 通道 中 被 RFI 损 坏 
的 原始 滤波 器 组 数据 . 在 给 定 的 时 间 样 本 中 , 识 
别 并 删除 比 中 值 多 3 倍 能 量 的 通道 . 使 用 已 知 色 
散 量 (Dispersion measure, DM) 56.8157 pe- cm-3, 
对 去 除 射 频 干 扰 后 的 数据 进行 非 相 干 消 色散 ,并 
将 其 压缩 到 一 个 频率 频道 . 大 部 分 脉冲 星 的 自转 
周期 是 非常 稳定 的 , 偶尔 会 在 短 时 间 内 加 速 , 这 
一 现象 被 称 为 周期 跃 变 (GlitchbjB0,，Crab 脉 冲 星 是 
一 颗 年 轻 的 脉冲 星 ,其 旋转 速率 会 发 生 Glitch[al. 
因此 , 利用 脉冲 星 计 时 包 TEMPO2 对 Crab 脉 冲 星 
到 达 NSRT 的 时 间 进 行 长 期 监测 9. 随后 , 利用 星 
历 表 确定 多 项 式 系 数 , 将 数据 按照 预测 的 周期 进 
行 单 脉冲 折合 . 在 整个 周期 内 , 脉冲 轮廓 由 1024 个 
相位 bin 构 成 . 对 脉冲 星 进 行 12.6 h 的 连续 观测 , 涵 
六 了 106 个 单 脉 冲 . 对 数据 进行 搜索 , 寻找 信 噪 比 
(Signal-to-Noise ratio, S/N) 阐 值 高 于 10 的 候选 体 . 
对 所 选 脉冲 的 子 集 进行 可 视 化 分 析 , 确认 所 有 窜 
带 RFI 都 被 移 除 . 

在 对 脉冲 星 观 测 的 12.6 h 中 , 共 探 测 到 2097 个 
信 噪 比 大 于 10 的 GP, 平均 每 21.63 s 探 测 到 1 个 GP， 
或 每 640 个 自转 周期 中 探测 到 1 个 GP. 如 上 文 所 述 ， 
GP 只 发 生 在 脉冲 星 旋转 的 MP 或 IP 相 位 内 回 . 图 1 
中 显示 了 两 个 GP 示例 , 其 中 脉冲 分 别 在 MP ( 左 ) 和 
IP ( 右 ) 辐 射 相位 中 被 检测 到 . 

从 图 1 可 以 看 出 , 频率 相位 图 像 中 未 显示 拖 尾 ， 
这 证 实 了 DM 的 值 是 正确 的 . 由 于 没有 进行 基于 品 
声 发 生 器 的 绝对 流量 测量 , 因此 利用 辐射 计 方程 估 
算 各 GP 的 流量 密度 B31, 方程 如 下 : 


(S/N) Ssys 
/Nptint AV 


FHS /NA TAR IL, Ssys = Ssyso + ScN 为 总 系统 
等 效 流量 密度 (System Equivalent Flux Density, 


S= (1) 


Tit 
SEFD), Ssyso = 200 Jy 是 天 空 背景 辐射 情况 下 的 


E ESE: Crab 脉 冲 星 巨 脉冲 辐射 等 待 时 间 分 布 的 非 稳 态 泊 松 模 拟 
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2 期 


SEFD, Son = 955v~°?? (vy 是 以 GHz 为 单位 的 频 
8) EHR ES ASEFDS4), np = 2 为 总 的 偏振 
数目 , tine = 32 Ls 为 采样 时 间 , Av = 342 MHz 为 
有 效 观 测 带宽 . 我 们 的 探测 阔 值 对 应 于 通 量 密度 为 
100 Jy. 
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上 对 间 序 列 图 


0.92 


Et 


图 1 分 别 在 MP ( 左 ) 和 IP ( 右 ) 辐 射 窗 
动态 谱 和 频率 全 加 后 的 


Fig.1 Dedispersed dynamic spectra and frequency averaged 
time series of two GPs detected in the MP (left) and IP 
(right) phases, respectively 


Rtas T PRIM BY AGP AB ic Bn Ze ak E 
电 辐 射流 量 的 直方 图 . 图 中 灰色 表示 Crab 脉 冲 星 
的 射电 积分 轮廓 . 该 图 清晰 表明 , 本 观测 所 获得 的 
GP 广 泛 分 布 于 积分 轮廓 所 定义 的 辐射 窗口 内 . 然 
而 , 前 人 的 多 项 观测 研究 发 现 Crab 脉 冲 星 的 GP 位 
于 积分 轮廓 的 相位 前 沿 B, 35-36], See Bera") Ail 
用 15 m 射 电 望 远 镜 在 1330 MHz 对 Crab 脉 冲 星 观 
测 了 约 260 h, 发 现 GP 的 轮廓 形状 具有 多 种 结 
构 特征 . 有 些 GP 有 单个 辐射 成 分 有些 GP 表 现 
出 多 个 复杂 的 辐射 成 分 . 而 且 GP 出 现 的 相位 
也 是 随机 的 , GP 位 于 积分 轮廓 的 相位 前 沿 、 中 
间 和 后 沿 没有 表现 出 任何 规律 性 . 从 图 2 看 
出 ， 共 有 1960 个 GP 位 于 辐射 的 MP 窗口 内 , 平均 
每 23.1 s 有 1 个 , 而 只 有 136 个 GP 位 于 IP 窗 口内 , 平 
均 每 332.5 s 探 测 到 1 个 . GP 搜 寻 通 常 采用 的 是 阔 
值 算法 , 即 寻 找 信 噪 比 高 于 某 个 阔 值 的 候选 体 . 
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Karuppusamy 等 站 利用 Westerbork 综 合 孔 径 射 电 
望远镜 对 Crab 脉 冲 星 的 GP 进行 观测 研究 , 采用 信 
We CG BY EL AFTRA GP, 共 探 测 到 12959 个 GP 爆 
发 ， 其 中 有 11384 个 GP 在 MP 相位 处 被 探测 到 ， 
1370 个 GP 在 IP 处 被 探测 到 . Berat BTA] H15 m 射 
电 望 远 镜 在 1330 MHz 对 Crab 脉 冲 星 的 GP 进行 了 
搜寻 , 采用 信 噪 比 阔 值 10, 共 探 测 到 1799 个 GP 爆发 
现象 , GP 的 能 量 服从 窜 律 谱 分 布 特征 , 谱 指 数 约 为 
一 3. 正常 脉冲 和 GP 辐射 具有 不 同 的 能 量 分 布 特征 ， 
GP 的 能 量 分 布 严格 遵从 究 律 分 布 B839], 然而 正常 
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探测 到 的 GP 在 MP ( 左 ) 和 IP ( 右 ) 窗 
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脉冲 的 能 量 分 布 服 从 对 数 正 态 分 布 B3; 和 9. 正常 肪 


冲 和 GP 能 量 分 布 不 同 可 能 是 由 它们 辐射 的 物理 
机 制 不 同 导致 的 . 为 了 验证 GP 搜寻 方法 , 保证 搜 
寻 结 果 不 会 对 后 续 的 科学 结论 产生 影响 , 我 们 计 
算 了 所 有 探测 到 的 GP 的 能 量 . 图 3 展示 的 是 探测 到 
的 2097 个 GP 能 量 统计 分 布 图 . GP HEB A AR 
律 谱 分 布 , 谱 指数 为 -3.5 + 0.2 (图 3 中 直线 ), 而 正 
常 脉冲 能 量 服从 对 数 正 态 分 布 , 二 者 有 本 质 上 的 区 
别 . 因此 采用 信 噪 比 阔 值 高 于 10 搜 寻 GP 的 方法 是 
可 行 的 , 不 会 对 后 续 的 科学 结论 产生 影响 . 
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Fig.2 Histograms showing the phase distributions of the detected GPs in the main pulse (left) and interpulse (right) emission 
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图 3 ”2097 个 GP 的 脉冲 能 量 统 计 分 布 
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Fig.3 Statistical distribution of the pulse energies for the 
2097 giant pulses 


3 结果 


3.1 “ 非 稳 态 的 泊 松 过 程 

研究 表明 , GP 的 辐射 机 制 是 在 具有 时 空 尺 度 
不 变性 的 自 组 织 临 界 状态 (selforganized criticality, 
SOC) FFER. 在 通常 情况 下 , 认为 GP 是 独立 
随机 发 生 的 事件 , 并 遵循 具有 指数 WTD 的 标准 泊 


松 过 程 : 


P(At) =e", (2) 


中 t 代 表 观 测 时 长 ,At 代表 连 续 GP 的 WT, AW 
GP 的 爆发 率 . 图 4 显示 了 5 min 内 GP 的 平均 爆发 率 
的 时 间 演 变 , 表明 GP 爆发 率 和 (t) 随 时 间 变 化 . GP 在 
MP 相位 和 IP 相 位 时 间 分 布 分 别 显 示 在 中 图 和 下 
K, GP 的 爆发 率 为 平均 每 5 min 一 次 . 平滑 滤波 后 
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的 分 布 用 蓝 色 实 心 曲线 表示 , 顶部 区 域 的 红 圈 标记 比 与 MP 辐射 相关 的 GP (GPs in the MP window, 

了 GP 生成 时 对 应 的 事件 时 间 . GMP) 更 为 显著 . GIP 的 辐射 在 低 波峰 率 ( 入 ) 的 时 间 

得 到 的 时 间 依 赖 或 非 平稳 泊 松 过 程 给 出 WTD 间隔 内 非常 活跃 , 根据 Wheatlandl43 的 工作 , 这 种 

为 聚集 现象 表明 GP 爆 发 率 J(A) 的 时 间 分 布 可 能 遵循 
JE MYPA DAt A 一 个 指数 函数 . 所 观察 到 的 分 布 可 近似 为 

So 和 Ddt F(A) = Are, (5) 


其 中 观测 时 间 间 隔 为 [0, T), 如 果 t > T, 则 出 现 率 
AT) = 0. 表达 式 简 化 为 
_ for f(A) Ne tdA 
EAFA ? 4) P(At) = ae- 2 (8 十 Ah“ 3， (6) 
EXS (AJAA = dt/T, 即 f( 和 ) 表 示 GP 爆 发 率 的 时 间 | 
分 布 . 在 观测 时 间 区 间 [0, TP, 平均 爆发 率 为 x 其 中 6 = [PQ 一 /入 ”和 0 < a < 2, Tyi 
= JOO). 时 变 泊 松 过 程 的 WTD (QRR A 根据 Guidorai 等 “的 工作 , (6) 式 可 以 化 简 为 
赖 于 GP 爆发 率 f( 和 ) 的 时 间 分 布 , GP 的 爆发 率 分 布 P(At) = (2 -a)8 (8+ At, (7) 
5K BEEBE A AM RAIA. 如 图 4 所 示 , 可 以 
清楚 地 看 到 以 瞬间 涨 落 为 特征 的 高 间 欣 性 , 这 反 其 中 a 决定 了 聚集 程度 , 6 表征 了 WTD 中 断 的 WT. 
映 了 GP 可 能 以 聚集 的 形式 出 现 . 特别 是 对 于 IP 阶 ”对 于 非常 大 的 WT (At > 6), CRE RIER IR, 
段 的 GP (GPs in the IP window, GIP), 聚集 效应 BAB — a. 


4 


时 中 A 是 归 一 化 项 , a 和 B 是 自由 参数 . 将 (5) 式 代 
ANARH, 得 到 WTD 的 解析 表达 式 为 


P(At) 


i ee fs se © mT 100. 


0 2 4 6 8 10 12 
Time/h 


图 4 探测 到 的 GP 在 整个 辐射 窗口 (上 )、MP 窗 口 ( 中 ) 和 JP 窗口 (下 ) 内 爆发 率 随时 间 的 演化 关系 


Fig.4 The temporal evolution of burst rate for detected GP in the whole emission window (upper), main pulse window (middle) 


and interpulse window (lower) 


19-5 
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3.2 GP 的 WTD 
将 Crab 脉 冲 星 探测 到 的 GP 的 WTD 如 图 5 所 

示 , 图 中 为 对 数 坐 标 . 对 巨 脉冲 的 WT 进行 统计 分 
ir, 横 坐 标的 间隔 为 100, 该 统计 分 布 中 WT 较 长 的 
数量 比较 少 . 应 用 非 平 稳 泊 松 分 布 的 模型 ((7) 式 ) 对 
其 进行 最 优 拟 合 , FUG He RB. 在 低 计 数 
的 情况 下 , 使 用 简单 的 卡 方 x 最 小 化 拟 合 每 个 bin 
中 收集 的 计数 , 其 结果 与 预期 分 布 不 符 , 取而代之 
的 是 可 以 采用 贝 叶 斯 马尔 可 夫 链 蒙特 卡 罗 (Markov 
Chain Monte Carlo) 算 法 . 值得 注意 的 是 , 对 于 大 
的 WT, 模型 与 观测 结果 之 间 存 在 偏差 , 这 可 能 是 1 
于 大 样本 具有 高 统计 灵敏 度 造成 的 . 对 于 长 WT 的 
差异 , 特别 是 对 于 低 爆 发 率 , 也 可 能 是 由 于 观测 确 
定 的 分 布 偏离 非 平 稳 泊 松 模型 . 
10° Em 一 一 一 一 一 
上 一 一 WTD of all GPs 了 
---: Non-stationary Poisson | 


10775 


Occurrence Probability 


1074} 
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45 来自 Crab 脉 冲 星 GP 的 WTD 图 和 利用 非 稳 态 泊 松 过 程 的 最 佳 

拟 合 结果 

Fig.5 The WTD of all GPs detected from the Crab pulsar 


with the corresponding best-fit model of a non-stationary 


Poisson process 


分 别 计算 连续 GMP 和 GIP 之 间 的 WT, GMP 
和 GIP 的 WTD 比 较 见 图 6. 在 GMP 和 GIP 之 间 可 以 
识别 出 非常 不 同 的 WTD 分 布 , 并 且 两 者 都 显示 出 
AMEE. 除了 GIP 的 窘 律 指数 明显 较 平坦 外 , GP 
和 GMP 均 可 以 进行 指数 拟 合 . 考虑 到 不 同类 型 的 
GP, 这 是 值得 注意 的 . 两 种 WTD 均 拟 合 为 非 平稳 
泊 松 分 布 ((7) 式 ). 表 1 给 出 了 所 有 WTD 的 最 佳 拟 合 
参数 . 估计 的 指数 a 和 6 差异 较 大 , 这 表明 IP 和 MP 


相 的 GP 可 能 来 自 不 同 的 辐射 机 制 |. 
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Waiting Time/s 
6 GP 分 别 在 MP 和 IP 辐 射 窗口 内 的 WTD 图 和 对 应 的 最 佳 拟 合 结 
果 
Fig.6 WTDs of GPs at the MP and IP phases with the 


corresponding best-fit models of a non-stationary Poisson 


PR 


process 
#1 利用 (7) 式 对 不 同 WTD 进 行 拟 合 的 最 佳 参数 
Table 1 Best-fit parameters obtained for 
different WTDs by using the model in Eq. (7) 


Sample Size a B 
GPs (all) 2096 0.207308 27.3974 8 
GPs (MP) 1960 0.18'§33 27.291151 


GPs (IP) 136 0.94+010 g4.45+1.98 


4 wit 

本 文 研究 了 Crab 脉 冲 星 12.6 hb 内 探测 到 GP 的 
WTD, 采用 非 平稳 泊 松 过 程 对 观测 数据 进行 了 贝 
叶 斯 统计 分 析 , 发 现 与 观测 到 的 WTD 是 一 致 的 , 结 
果 表 明 GP 是 独立 的 随机 事件 . GP 的 爆发 紊 在 一 定 
程度 上 反映 了 能 量 的 释放 速率 . 一 般 认 为 脉冲 星 主 
要 以 磁 侦 极 辐 射 的 形式 损失 能 量 , 这 导致 其 旋转 速 
度 相应 减 慢 , 这 可 能 是 由 于 固定 涡 度 突变 释放 产生 
pua, 由 于 超 流 体外 核 和 地 过 之 间 的 角 动量 交换 
可 以 改变 脉冲 星 的 磁 层 . 因此 GP 也 可 能 与 脉冲 星 
磁 层 中 的 随机 动力 学 过 程 有 关 . Palfreyman 等 [和 9 发 
现 , 在 船 帆 座 脉冲 星 发 生 小 的 周期 跃 变 后 , GP 爆发 
率 也 相应 增加 . 此 外 , 在 PSR J0742-2822 中 显示 辐 
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射 的 模式 变换 与 Glitch 之 间 也 存在 联系 6. Crab 肪 


系 . 


盾 


冲 星 的 Glitch 是 非常 活跃 的 的 ,进行 长 期 监测 非常 
有 必要 , 这 样 就 可 以 揭示 GP 的 爆发 率 与 Glitch 的 联 


观测 结果 确定 WTD 与 SOCI 的 解释 并 不 矛 
并 且 与 GP 的 非 平稳 模型 一 致 . 当 爆 发 率 变化 


时 , 可 以 再 现 观测 到 的 WTD 的 定性 特征 . 另 一 种 解 


释 是 间歇 性 等 离子 体 消 流 驱动 Crab 脉 冲 星 的 GP. 
然而 到 目前 为 止 , 还 没有 理论 模型 能 够 正确 预测 


WT 


D. 
本 研究 分 别 考虑 了 GMP 和 GIP 的 WT, 并 给 出 


了 不 同 的 分 布 . GIP 比 GMP 更 具有 聚集 性 . 一 般 来 
W, MP 和 IP 辐 射 被 认为 是 来 自 星体 的 两 个 磁极 . 如 


果 辐 射 区 域 的 物理 条 件 相似 , 那么 在 MP 和 IP 中 活 


跃 的 射电 辐射 机 制 就 是 相同 的 . 观测 确定 的 GMP 


和 G 
EG 


IP 的 WTD 有 显著 差异 , 这 表明 在 相反 的 磁极 
P 的 辐射 机 制 是 不 同 的 . 


此 外 , 通过 对 Crab 脉 冲 星 [9 中 个 别 GP 的 高 时 


间 分 辩 率 观测 , 发 现 GIP 的 特征 在 时 域 和 频 域 上 都 
与 GMP 的 特征 有 明显 不 同 . GMP 偶 尔 可 以 分 解 为 
短暂 的 和 相对 罕 带 的 爆发 , 这 一 现象 是 由 强 等 离子 


体 满 流 中 的 孤子 击 缩 产生 的 . 射电 发 射 区 的 等 离子 
体 流动 是 动态 的 . GIP 的 特征 是 相对 狭窄 、 微 秒 级 


的 发 射 带 . 这 些 按 比 例 间隔 的 波段 几乎 同时 开始 ， 
偶尔 在 频率 上 有 轻微 的 上 升 . 因为 每 个 模型 都 预测 


了 在 等 离子 体 频 率 下 有 规则 间隔 的 窗 带 辐射 , 所 
以 己 知 的 脉冲 星 射电 辐射 模型 难以 解释 IP 的 高 动 
态 频 谱 . 参考 太阳 泡 斑 的 物理 特性 , 可 以 应 用 共振 
回旋 辐射 和 几何 效应 这 两 类 模型 来 解释 GIP 的 剧 


烈 动 态 谱 区 .根据 共振 回旋 辐射 模型 外 ], 在 脉冲 星 
高 纬度 区 域 , GIP 是 由 靠近 光 柱 的 高 能 粒子 在 磁场 
中 受 力 加 速 产 生 的 . 对 几何 效应 模型 , 辐射 光束 被 
相干 地 分 裂 并 与 自身 干涉 , 从 而 产生 干涉 条 纹 [9]. 
GIP 和 GMP 之 间 的 显著 差异 表明 它们 在 脉冲 星 磁 
层 内 的 发 射 机 制 或 传播 效应 可 能 不 同 . 由 于 GP 在 


更 高 频率 的 IP 相 位 更 加 常见 , 所 以 为 了 揭示 GP 的 


物理 


辐射 机 制 , 需要 在 更 高 频率 、 更 大 带宽 和 更 短 


时 间 分 辨 率 条 件 下 观测 . 


WT 的 统计 性 质 对 于 理解 脉冲 星 磁 层 的 物理 性 


质 非常 有 用 , 在 少数 具有 大 磁场 值 的 高 能 脉冲 星 中 


CAG Sl SGP HOO ON, 对 它们 进行 长 期 连续 的 


监测 , 可 以 发 现 它 们 共有 相同 或 不 同 的 分 布 . 


5 结论 

本 文 研究 了 Crab 脉 冲 星 GP 的 WT 统计 . 我 们 
对 Crab 脉 冲 星 进 行 了 长 时 间 连 续 观 测 ，GMP 和 
GIP 的 WTD 都 具有 聚集 特征 , 且 符 合 非 平 稳 泊 松 
过 程 , 其 中 GIP 的 爆发 率 具 有 高 度 波 动 性 . GMP 和 
GIP 的 WTD 之 间 存 在 明显 差异 , 这 表明 它们 可 能 
源 于 不 同 的 物理 辐射 机 制 . 
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Apstract We study the waiting time distribution of giant pulses from the Crab pulsar observed with 
the Nanshan 26-m radio telescope at a center frequency of 1556 MHz with 512 MHz bandwidth. The 
observations were performed over a duration of 12.6 hours with 32 us sampling interval. Our analysis 
has led to the detection of 2097 giant pulses above a threshold of 10 o, with flux density > 100 Jy. 
The occurrence rate of giant pulses is characterized by a highly intermittent giant pulse productivity in 
short clusters with high rates, separated by relatively long quiescent intervals with low occurrence rates, 
especially for the giant pulses associated with the interpulse emission. The distribution of waiting times 
between two subsequent giant pulses displays a power-law tail and can be modeled with a non-stationary 
Poisson process, which indicates that giant pulses are independent and random events. Distinct waiting 
time distributions between giant pulses in the main pulse and interpulse phases are presented, which 
implies that the giant pulse emission mechanisms maybe different in the opposite magnetic poles. The 
ramification for our understanding of the radio emission mechanisms is discussed. 
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